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ABSTRACT 

We present Chandra studies of the X-ray binary (XRB) populations in the bulge and ring regions of the 
ring galaxy NGC 1291. We detect 169 X-ray point sources in the galaxy, 75 in the bulge and 71 in the 
ring, utilizing the four available Chandra observations totaling an effective exposure of 1 79 ks. We report 
photometric properties of these sources in a point-source catalog. There are « 40% of the bulge sources and 
% 25% of the ring sources showing > 3a long-term variability in their X-ray count rate. The X-ray colors 
suggest that a significant fraction of the bulge (« 75%) and ring (« 65%) sources are likely low-mass X-ray 
binaries (LMXBs). The spectra of the nuclear source indicate that it is a low-luminosity AGN with moderate 
obscuration; spectral variability is observed between individual observations. We construct 0. 3-8.0 keV X- 
ray luminosity functions (XLFs) for the bulge and ring XRB populations, taking into account the detection 
incompleteness and background AGN contamination. We reach 90% completeness limits of ~ 1.5 x 10 37 and 
~ 2.2 x 10 37 erg s~ l for the bulge and ring populations, respectively. Both XLFs can be fit with a broken 
power-law model, and the shapes are consistent with those expected for populations dominated by LMXBs. 

We perform detailed population synthesis modeling of the XRB populations in NGC 1291, which suggests 
that the observed combined XLF is dominated by an old LMXB population. We compare the bulge and ring 
XRB populations, and argue that the ring XRBs are associated with a younger stellar population than the bulge 
sources, based on the relative overdensity of X-ray sources in the ring, the generally harder X-ray color of the 
ring sources, the overabundance of luminous sources in the combined XLF, and the flatter shape of the ring 
XLF. 

Subject headings: galaxies: individual (NGC 1291) — galaxies: luminosity function — X-rays: binaries 
X-rays: galaxies 


1. INTRODUCTION 

Galactic rings are a peculiar phenomenon among galaxies. 
A prototype ring galaxy is the famous Cartwheel galaxy, with 
a spectacular ring luminous in the infrared (IR), optical, UV, 
and X-ray bands. The formation of ring-like structures in 
galaxies is not fully understood. The two most common theo- 
ries are the collisional origin and the resonance origin. A col- 
lisional ring could be created by a head-on collision between 
a compact companion galaxy and a larger disk system; a col- 
lisional ring galaxy usually contains an off-centered nucleus 
and an asymmetric ring structure, and there is a significant 
chance of finding a companion galaxy close to the ring galaxy 
(see Appleton & Struck-Marcell 1996 for a review). A reso- 
nance ring could be produced by the outer Lindblad resonance 
of spiral bars; such a system is preferentially observed in 
barred galaxies, and the resulting ring is generally symmetric 
(see Buia & Combes 1996 for a review). There are a few sys- 
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terns that are difficult to explain with these two theories, such 
as the Hoag’s Object (Hoag 1950), for which alternative for- 
mation scenarios have been proposed (e.g., Finkelman et al. 
201 1). 

Galactic rings are expected to exhibit significant star for- 
mation compared to their progenitor galaxies mainly due to 
the enhanced gas densities in the rings, regardless of the 
formation mechanism. Surveys of ring galaxies in optical, 
Ha, IR, and radio have revealed that the rings arc generally 
rich in Hi gas and Hll regions (e.g., van Driel et al. 1988; 
Marston & Appleton 1995; Phillips 1996; Ryder etal. 1996; 
Grouchy etal. 2010). With Chandra's subarcsecond angular 
resolution, it is now possible to study in detail the X-ray bi- 
nary (XRB) populations in nearby galaxies, e.g., low-mass 
X-ray binaries (LMXBs) that are associated with old stellar 
populations and high-mass X-ray binaries (HMXBs) that are 
associated with young stellar populations (see e.g., Fabbiano 
2006 for a review). Therefore, X-ray observations have be- 
come an independent approach for constraining the nature of 
the underlying stellar populations by means of X-ray lumi- 
nosity functions (XLFs) and X-ray colors of the XRBs. How- 
ever, little is known regarding the X-ray properties of galac- 
tic rings, except for the Cartwheel galaxy, which has an ex- 
ceptionally large number of young (< 10 7 yr) ultraluminous 
X-ray sources (ULXs; L x > 10 39 erg s -1 ; e.g., Swartz etal. 
2004) in the star-forming ring (e.g., Gao et al. 2003; King 
2004; Wolter & Trinchieri 2004). 

NGC 1291 is a remarkable ring galaxy, first noted 
by Perrine (1922). It was classified as (R)SB(s)0/a 
(RC3; de Vaucouleurs et al. 1991), and is virtually face-on 
(Meboldetal. 1979). The images of NGC 1291 in the UV, 
optical, and near-IR (NIR) are shown in Figure 1. The 
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galaxy includes a fainter outer ring and a prominent bulge, 
which was considered as a prototypical pseudobulge (e.g., 
Kormendy & Kennicutt 2004). At a distance of 8.9 Mpc, 
NGC 1291 has Af B) o— 20.6 (de Vaucouleurs 1975). The D 25 
ellipse of NGC 1291 has a size of 9.8' x 8.1' (Kennicutt et al. 
2003), and the nuclear position of the galaxy is aj2ooo o = 
03 h 17 m 18.6 8 , <Sj 2 ooo.o = — 41°06'29" (Evans et al. 2010). The 
Galactic column density along the line of sight to NGC 1291 
is Nr =2.1 x 10 20 cm"* (Dickey & Lockman 1990). 10 

The ring of NGC 1291 does not appear as energetic as that 
of the Cartwheel galaxy. There were only two ULXs re- 
ported in the ring based on the first Chandra observation in 
2000 (Swartz et al. 2004). The total star formation rate (SFR) 
of NGC 1291 is only 0.4 A/ 0 yr~ L (Kennicutt et al. 2003), 
compared to that of 67.5 Mq yr~ ] for the Cartwheel galaxy 
(Higdon 1993). Given its barred spiral nature, the relatively 
symmetric structure of the ring, and the absence of any ap- 
parent companion, NGC 1291 is likely a resonance ring sys- 
tem, while the Cartwheel galaxy is most likely a collisional 
ring system (e.g., Appleton & Struck-Marcell 1996). Intense 
star formation still appears present in the ring of NGC 1291, 
which is obvious from the UV image (Figure 1). On the con- 
trary, the bulge is much redder than the ring, and dominates 
in the NTR image (Figure 1). van Driel et al. (1988) mapped 
the Hi distribution in NGC 1291, and revealed that Hi gas is 
concentrated in the ring region. The ring also contains dis- 
crete Hu regions, while the bulge has filaments of ionized 
gas (Crocker et al. 1996). Given the UV/optical/radio data, 
the ring appears to host a younger stellar population than the 
bulge. Noll et al. (2009) studied the star-formation history of 
many galaxies including NGC 1291, using spectral energy 
distribution fitting techniques. The best-fit model for NGC 
1291 (Noll, private communication) includes two exponential 
starburst events: one happened 10 Gyr ago with an e-folding 
timescale of 0.25 Gyr and the other happened 200 Myr ago 
with an e-folding timescale of 50 Myr. The stellar mass of the 
old population is « 330 times higher than that of the younger 
one. Although the fitting results were derived for the entire 
galaxy and may have large uncertainties, they are consistent 
with the UV/optical/radio data if the ring hosts the younger 
stellar population. 

NGC 1291 had been observed with Chandra twice in June 
and November 2000, with a total exposure of « 80 ks. A first 
study by Irwin et al. (2002) was limited to the bulge region. 
Kilgardet al. (2005) presented a point-source catalog based 
on these observations, and found that the XLF and X-ray col- 
ors suggest old XRB populations in the galaxy. This study 
however, did not discriminate between bulge and ring XRBs. 

We obtained two additional Chandra observations of NGC 
1 29 1 in May 2010, with a total exposure of « 1 20 ks. The new 
observations cover the majority of the ring region, allowing 
us to investigate the XRB populations in the bulge and ring 
respectively. NGC 1291 is much closer than the Cartwheel 
galaxy (« 10 Mpc compared to « 120 Mpc), and thus we are 
able to probe the ring XRB population down to an unprece- 
dented luminosity limit (« 1 0 37 erg s -1 ). We compare the ring 
XRB population to the bulge population and test whether the 
galactic ring hosts a younger stellar population from an X-ray 
point of view. The deep observations on the bulge also allow 
us to compare its LMXB population to those of the recently 
acquired deep exposures of elliptical galaxies, and establish if 
all LMXB populations are indeed compatible. 

10 See http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl. 


This paper is organized as follows. We describe the Chan- 
dra observations of NGC 1291 and the data reduction proce- 
dure in Section 2. In Section 3 we construct an X-ray point- 
source catalog. X-ray source properties are presented in Sec- 
tion 4, including X-ray source radial profile, source variabil- 
ity, and X-ray colors. The nuclear source and ULX candidates 
are also discussed. In Section 5, we calculate and parameter- 
ize the XLFs for the bulge and ring regions, respectively, tak- 
ing into account the detection incompleteness and background 
AGN contamination. We model the observed XLF using pop- 
ulation synthesis (PS) models in Section 6. In Section 7, we 
discuss implications of the X-ray properties and we argue that 
the ring hosts a younger stellar population than the bulge. We 
summarize our results in Section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 1291 has been observed with the S3 chip of Chandra 
Advanced CCD Imaging Spectrometer (ACTS; Garmire et al. 
2003) at four epochs. The four observations are listed in Ta- 
ble 1 along with their exposure times, ranging from 37 ks to 
70 ks. The focal-plane temperature was - 1 20°C for all the ob- 
servations. We reduced and analyzed the observational data 
using mainly the Chandra Interactive Analysis of Observa- 
tions (CLAO) tools. 11 We used the chandrarepro script 
to reprocess the data with the latest calibration. The back- 
ground light curve of each observation was then inspected and 
background flares were removed using the DEFLARE CIAO 
script. The flare-cleaned exposure times are also listed in Ta- 
ble 1; there was a significant flare in observation 2059 that 
reduces its usable exposure by ~ 40%. The total usable expo- 
sure for NGC 1291 is 178.6 ks. 

We registered the astrometry of the observations using the 
CIAO tool REPROJECT_ASPECT. We created a 0.3-8.0 keV 
image for each observation and searched for sources using 
wavdetect (Freeman et al. 2002) at a false-positive proba- 
bility threshold of 10~ 6 . These source lists were used to reg- 
ister the observations to the astrometric frame of observation 
11272 with the CIAO script REPROJECT ASPECT, adopting 
a 3" matching radius and a residual rejection limit of 0.6". 
We reprojected the registered event lists to observation 1 1272 
using REPROJECT_EVENTS, and then merged all the obser- 
vations to create a master event file using dmmerge. The 
ACIS-S3 chip has different pointings for the four observa- 
tions, and the average aim point (weighted by exposure time) 
is a, 2000.0 = 03 h 17 ra 09.31 s , <Wo = -41°06 , 25.1". 

We constructed images from the merged event file using 
the standard ASCA grade set ( ASCA grades 0, 2, 3, 4, 6) 
for six bands: 0.3-8.0kcV (fiill band), 0.3-2.0keV (soft 
band), 0. 5-2.0 keV (conventional soft band), 2.0-8.0keV 
(hard band), 0.3-1 .0 keV (soft band 1), and 1 .0-2.0 keV (soft 
band 2). For each observation, we created exposure maps in 
these bands following the basic procedure outlined in Sec- 
tion 3.2 of Hornschemeier et al. (2001), which takes into ac- 
count the effects of vignetting, gaps between the CCDs, bad- 
column filtering, bad-pixel filtering, and the spatially depen- 
dent degradation in quantum efficiency due to contamina- 
tion on the ACIS optical-blocking filters. A photon index of 
F = 1.7 was assumed in creating the exposure map, which 
is a typical value for photon indices of X-ray binaries (e.g., 
Irwin et al. 2003; Brassington et al. 2010). Merged exposure 
maps were then created from the exposure maps of the indi- 
vidual observations. 

1 1 See http://cxc.harvard.edu/ciao/ for details on CIAO. 
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Fig. 1. — Images of NGC 1291 in the GALEX ncar-UV (NUV, ^ = 2267 A), Digitized Sky Survey red (RG610), and 2MASS AT bands. The size of each 
image is 14' x 14', centered on the galactic nucleus. The ellipses outline the bulge and ring regions (see Section 3 for details). The UV ring appears to be 
dominated by young stars. 
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Right Ascension (2000) 

Fig. 2. — Chandra “false-color” image of NGC 129 1 . This image is a color 
composite of the exposure-corrected adaptively smoothed images in the 0.3- 
1 .0 keV (red), 1. 0-2.0 keV (green), and 2.0-8.0 keV (blue) bands. The cross 
symbol shows the center of the galaxy. The white dashed ellipse shows the 
D 25 region (Kennicutt et al. 2003). The cyan and magenta ellipses outline the 
bulge anc ring regions of the galaxy, respectively. 


We created adaptively smoothed images from the raw im- 
ages using the Cl AO tool CSMOOTH. Exposure-corrected 
smoothed images were then constructed following Section 3.3 
of Baganoff et al. (2003). We show in Figure 2 a color com- 
posite of the exposure-corrected smoothed images in the 0.3- 
1.0, 1. 0-2.0, and 2.0-8.0 keV bands. 

3. POINT-SOURCE CATALOG 


We ran WAVDETECT on the full-band image 12 to search for 
X-ray point sources, with a “\/2 sequence” of wavelet scales 
(i.e., 1, \/2, 2, 2y/2, 4, 4\/2, 8, 8\/2, and 16 pixels) and a 
false-positive probability threshold of 1 0 -6 . At this threshold, 
we expect approximately two false detections given the size 
of the image (e.g., Kim et al. 2004). 

There are 169 X-ray sources detected in the full band. 
We utilized the ACIS EXTRACT (AE; Broos et al. 2010) pro- 
gram to derive the source photometry. The main advantage 
of AE compared to a simple circular-aperture extraction ap- 
proach is that AE uses a polygonal extraction region to rep- 
resent the point spread function (PSF) for each observation 
of each source individually, which is more suitable for mul- 
tiple observations with different aim points and roll angles. 
Briefly, AE constructs source extraction regions that approxi- 
mate » 90% encircled-energy fraction (EEF) contours of the 
PSFs at 1 .50 keV for each exposure. Smaller extraction re- 
gions (as low as ^ 60%) were used for sources in crowded 
regions. Background counts were extracted in an annular re- 
gion around each source excluding any pixels that belong to 
detected sources; the number of extracted background counts 
was required to be at least 100 to reduce the uncertainty from 
the background estimation. Encircled-energy fractions at sev- 
eral energies (ranging from 0.28 to 8.60 keV) were provided 
by AE, and thus we derived aperture corrections for ail the six 
X-ray bands listed in Section 2 via interpolation. Aperture- 
corrected net (background-subtracted) source counts and the 
corresponding ler errors (Gehrels 1986) for these bands were 
then computed; the errors were propagated through the errors 
of the extracted sources and background counts following the 
numerical method described in §1.7.3 of Lyons (1991). As 
the source detection was performed only in the full band, we 
consider a source as undetected in a given band (other than the 
full band) if it is less significant than 3cr, when calculated us- 
ing the net counts and la error in that band. An upper limit on 
the source counts was derived for an undetected source, using 
the Bayesian approach of Kraft et al. (1991) for the 99% con- 

12 Source searching was performed only in the full band as we focused on 
the study of the full-band XLFs. 
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fidence level if the extracted number of counts is less than 10 
in the given band; otherwise, we calculated a 3a upper limit 
following the Poisson statistic (Gehrels 1986). 

An X-ray spectrum for each source was extracted by AE 
along with the spectrum of its local background. We per- 
formed spectral fitting of each source using XSPEC (Version 
12.6.0; Arnaud 1996) with the Cash fitting statistic (CSTAT; 
Cash 1979). An absorbed power-law model (TBABS*POW) 
was employed for sources with more than 30 full-band counts; 
for the other sources, the spectra were modeled with a simple 
power law with the photon index fixed at F = 1 .7. The ob- 
served full-band fluxes were derived from the best-fit spectra, 
which were then converted to full-band luminosities. We ap- 
plied a smaller correction factor (1.05) to the observed full- 
band luminosities to account for the Galactic absorption. 

Of the 169 point sources, 75 are in the bulge region and 
71 in the ring area. We defined the bulge and ring regions 
based on the UV and NIR images shown in Figure 1. We 
examined the surface brightness contours of the A'-band im- 
age and selected a contour that encloses ~ 80% of the light 
within the D25 ellipse. This contour can be well approxi- 
mated by an 4.7' x 3.7' ellipse centered on the galactic nu- 
cleus, which was used to defined the bulge boundary. The ring 
region is defined by two ellipses centered on the galactic nu- 
cleus. We fixed the position angle (70°) and axis ratio (1.15) 
of the ellipses according to the morphology information given 
in de Vaucouleurs (1975) and only allowed the sizes of the el- 
lipses to vary so that the ring region includes the majority of 
the prominent UV ring features. We note that « 85% of ring 
region is covered by the Chandra exposure (see Figure 2). 

The spatial distribution of the X-ray sources is displayed in 
Figure 3a, and their luminosity distribution is show in Fig- 
ure 3b. The full-band luminosities range from 2.7 x 10 36 to 
2.3 x 10 39 erg s~ ! . We listed basic photometric properties for 
the 169 sources in Table 2, with details of the columns given 
below. 

1 . Column 1 : the source number, listed in order of increas- 
ing right ascension. 

2. Columns 2 and 3: the right ascension and declination 
of the X-ray source, respectively. 

3. Column 4: the off-axis angle of the X-ray source in ar- 
cminutes. This is the distance between the source posi- 
tion given in columns 2 and 3 and the average aim point 
(see Section 2). 

4. Column 5: the radial distance of the X-ray source to the 
nucleus of the galaxy in areminutes. 

5. Columns 6-23: the source counts and the corre- 
sponding la statistical errors or the upper limits 
on source counts for the six bands, 0.3-8.0keV, 
0.3-2.0 keV, 0.5-2.0keV, 2.0-8.0keV, 0.3-1. OkeV, 
and 1.0-2. OkeV, respectively. Upper limits are indi- 
cated as a “-1 .0” in the error columns. 

6. Column 24: the full-band effective exposure time de- 
termined from the exposure map, which has been cor- 
rected for vignetting and quantum-efficiency degrada- 
tion. 

7. Columns 25-30: the soft and hard X-ray colors and the 
corresponding la statistical errors; see Section 4.4 be- 
low for details. 


8. Column 31: the photon index derived from spectral fit- 
ting, or set to T = 1 .7 if the number of full-band counts 
is smaller than 30. 

9. Column 32: observed full-band X-ray luminosity, in 
units of erg s” 1 . 

10. Column 33: the statistical significance of the full-band 
flux variability; see Section 4.3 below for details. It is 
set to 1 -1 .0” if the flux variability is not available. 

11. Column 34: region flag for the source. The flag is 
“Bulge” if the source is located in the bulge, “Ring” 
if the source is in the ring, or “Other” if the source is 
neither in the bulge nor in the ring. 


4. POINT-SOURCE PROPERTIES 
4.1. Radial Profile 

We constructed a radial profile of the X-ray source surface 
density. Given the non-circular shapes of the ring and bulge, 
we calculated X-ray source densities in four elliptical annular 
regions in the bulge, four elliptical annular regions in the ring, 
and one region in between the bulge and the ring. The ellip- 
tical annular regions in the ring/bulge have the same position 
angle and axis ratio as the ring/bulge and are evenly spaced 
in their , distances to the galactic nucleus. We corrected the 
number of X-ray sources in each region for the detection in- 
completeness and background AGN contamination (see Sec- 
tion 5.1 below; the incompleteness correction was derived in 
the corresponding region). For this purpose, we only included 
sources that are outside the central 1 0"-radius region and have 
full-band luminosity L\ > 1 x 10 37 erg s -1 . The corrected 
number of sources was then divided by the region area to ob- 
tain the surface density. The X-ray source radial profile is 
presented in Figure 4. There appears to be an overdensity 
of X-ray sources in the ring of the galaxy compared to the 
region in between the bulge and the ring. We also included 
in this plot surface brightness radial profiles for the GALEX 
NUV and 2MASS £-band data, calculated in the same re- 
gions as above. The surface brightness profiles are arbitrarily 
normalized so that the highest brightness has the same value 
as the highest source density. The NUV surface brightness 
is expected to trace the young stellar populations, while the 
AT-band emission is sensitive to the old stellar populations. In 
the ring region, the NUV and X-band surface brightness pro- 
files differ significantly. The X-ray source radial profile in the 
ring does not continue with the trend shown in the bulge and 
the region in between, and it mainly follows the NUV surface 
brightness distribution, suggesting that the X-ray sources in 
the bulge and ring regions are associated with different stellar 
populations and the ring stellar population is likely younger. 

4.2. Central Point Source and XJltr aluminous X-Ray Sources 

The nuclear point source (XID 84) of NGC 1291 has a po- 
sition of aj2ooo.o - 03 h 17 ra l 8.58 s , feooo.o = -41 O 06'28.7", and 
an observed fhll-band luminosity of 2.3 x 10 39 erg s~ l . It has 
been classified as a low-ionization nuclear emission-line re- 
gion (LINER; e.g., Smith et al. 2007), and Irwin et al. (2002) 
suggested that it is an obscured low-luminosity active galac- 
tic nucleus (AGN) based on the X-ray spectral analysis of the 
first two Chandra observations. We show the X-ray spectra 
of the source for the four individual observations in Figure 5, 
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boundary lines show's the ACIS-S3 coverage of NGC 1291. (b) Full-band luminosity vs. distance to the nucleus of the galaxy. 



Fig. 4. — Radial profile of the X-ray source surface density. The densi- 
ties were computed in four elliptical annular regions in the bulge (green filled 
circles), four elliptical annular regions in the ring (magenta open circles), 
and one region in between the bulge and the ring (black triangle). The x-axis 
shows the average distance of each region to the galactic nucleus. The surface 
densities have been corrected for detection incompleteness and background 
AGN contamination. Only sources outside the central 10" -radius region and 
having L\ > 10 37 erg s” 1 are included in the density calculation. The errors 
of the surface densities were calculated based on the Poisson errors of the 
number of sources in each bin (Gehrcls 1986). The red and blue curves dis- 
play the surface brightness radial profiles (scaled to the highest X-ray source 
density) for the GALEXtWSV and 2MASS /f-band data, respectively. There 
appears to be an overdensity of X-ray sources in the ring region. (A color 
version of this figure is available in the online journal.) 

which were extracted by AE with polygonal regions that ap- 
proximate the « 89% EEF PSFs in the 0.3-8. 0 keV band. Lo- 
cal background spectra were also extracted by AE in the same 
annular regions as those used in the photometry extraction. 
There are clearly strong soft X-ray excesses in the spectra of 
all the observations, which are typical among low-luminosity 
AGNs and are considered to originate from hot gas in the 


galactic center (e.g., Ptak et al. 1999). We thus fit the spectra 
using XSPEC with an absorbed power-law (AGN) plus ther- 
mal plasma (hot gas) model (WABSl *APEC+WABS2*POW). 
Note that there could be multiple temperatures present in the 
gas, however, the limited signal-to-noise ratio (S/N) of the 
data cannot place strong constraints on a multi-temperature 
plasma. The spectral slope and column density for the power- 
law component were free to vary during the fitting. For the 
thermal component, we set the abundance to be Solar, and 
then we determined its temperature (0.15^;^ keV) and ab- 
sorption column density (Ah,i =7.6^7 x 10 21 cm -2 ) by fitting 
the stacked spectrum of the four observations. In the spectral 
fitting of the individual observations, the temperature and col- 
umn density parameters were fixed, and only the normaliza- 
tion was allowed to vary. 13 The best-fit parameters are listed 
in Table 3; also shown are the observed luminosities of the 
AGN component. The Id errors were derived by varying one 
interesting parameter. The intrinsic 2. 0-8.0 keV X-ray lumi- 
nosity is « (2-3) x 10 39 erg s -1 after absorption and aperture 
corrections, indicating its low-luminosity nature. The nuclear 
AGN has an intrinsic power-law photon index of ~ 2 and an 
absorption column density of « 2 x 10 22 cm -2 , which are typ- 
ical for an obscured AGN (e.g., Turner et al. 1997). The 3d 
confidence contours for the photon index and column density 
parameters are shown in Figure 6; there is a « 3a spectral 
variation between observations 795 and 1 1272. The X-ray lu- 
minosities of these two observations also differ at the « 3d 
significance level (Table 1). It has been noted that long-term 
variability in luminosity and spectral parameters is common 
among obscured AGNs (e.g., Risaliti et al. 2002). 

Besides the nuclear source, there are two more sources 
(XIDs 55 and 57) with observed full-band luminosity in ex- 

13 The extraction regions are different between observations, correspond- 
ing to the « 89% EEF PSFs of the nuclear point source in individual obser- 
vations. Therefore, the amount of diffused gas within the extraction region 
varies. 
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Fig. 5 . — X-ray spectra of the nuclear source overlaid with the best-fit mod- 
els for the four individual observations. The spectra were fit with an absorbed 
power iaw plus thermal plasma model (WABSl*APECH wabs2*pow). The 
dashed curves show the individual model components. For the thermal com- 
ponent, the temperature (0.15 keV) and column density (7.6 x 10 21 cm -2 ) 
were fixed at the values deri\ ed from the stacked spectrum. See Section 4.2 
for more details and see Table 3 for the fit parameters 

cess of 10 39 erg s’ 1 , in the regime of ULXs. Their in- 
trinsic full-band luminosities after absorption corrections are 
2.1 x 10 39 erg s" 1 and 2.0 x 10 39 erg s , respectively. One 
other source (XID 26), after being corrected for intrinsic ab- 
sorption (Nh = 1 .4 x 10 21 cm -2 ), has an intrinsic full-band lu- 
minosity of 1.1 x 10 39 erg s’ 1 , which should also be consid- 
ered as a ULX candidate. XIDs 26 and 55 are ring sources, 
whereas XID 57 is a bulge source. The two ring ULX candi- 
dates have been reported in the ULX catalog of Swartz et al. 
(2004) based on the data of the first observation (observa- 
tion 795). The bulge ULX candidate is a transient source, 
not detected in the first two observations. These ULX can- 
didates were not detected in the study of Irwin et al. (2002), 
which utilized only the first observation and focused on the 
bulge region. There is also a non-negiigible possibility that 
these sources are background AGNs. The expected numbers 
of background AGNs with ULX luminosities if placed at the 
distance of NGC 1291 are ~ 0.05 and ~ 0.3 for the bulge and 
ring regions, respectively (see Section 5.1 below for details 
of the AGN surface density). The numbers of ULXs in NGC 
1 291 and the Cartwheel galaxy (« 3 and ^ 20) do not appear 
to be linearly correlated with their SFRs (0.4 Af© yr’* and 
67.5 M 0 yr -1 ). Previous ULX surveys have found that the 
numbers of ULXs at low SFRs (< 0.4 M© yr“ l ) are signif- 
icantly higher than those extrapolated linearly from the high 
SFR end, suggestive of an additional population of ULXs as- 
sociated with LMXBs (e.g., Liu et al. 2006). 

4.3. Source Variability 

X-ray fiux/spectral variability is a common feature among 
point sources in galaxies, which is generally attributed 
to the change of physical properties of the accretion 
disks (e.g., Remillard & McClintock 2006; Done et al. 2007; 
Brassington et al. 2010; Fabbiano et al. 2010). The four 
Chandra observations of NGC 1291 span « 10 years, al- 
lowing us to study the long-term flux variability of the point 
sources. 

For each source in each observation, we calculated its 
aperture-corrected full-band count rate ( CR ) based on the pho- 
tometry provided by AE; an upper limit on the count rate was 


Fig. 6. — 3 a confidence contours for the photon index and absorption 
column density parameters of the nuclear AGN in the three individual ob- 
servations, observations 795 (blue), 1 1272 (green), and 1 1273 (red). The fit 
parameters in observation 2059 have large uncertainties and is not shown in 
the plot for clarity. The crosses show the best-fit A/h ,2 and T values in obser- 
vations 795, 1 1272, and 1 1273, respectively. The larger confidence contour 
of observation 1 1273 than that of observation 795 despite of its longer expo- 
sure is due to the lower fiux/luminosity of the nucleus during this observation. 
See Section 4.2 for more details. (A color version of this figure is available 
in the online journal.) 

calculated if the source was not detected in the observation 
(see Section 3). For a source that was detected in at least one 
individual observation and was covered by at least two obser- 
vations, we computed the maximum statistical significance of 
its flux variation between any two observations, defined as 
(e.g., Brassington et al. 2009; Sell et al. 201 1): 

\CR,-CRj\ 

o var = max, v ■ , ( 1 ) 

\J a CR, + "cRj 

where the subscripts i and j run over different observations, 
and o C R, and o C Rj are the ler errors of the count rates. We con- 
sider a source to oe variable if its o var parameter is greater than 
three (i.e., > 3 o variation). The flux variability for the point 
sources in NGC 1291 are shown in Figure 7. In the bulge 
and ring regions, 38^ lo /o and 25^ 3 % of the sources are vari- 
able, respectively; the errors are the la Poisson uncertainties 
on the numbers of sources that are considered variable. We 
note that the central point source and the three ULX candi- 
dates are variable. There is no significant difference between 
the fraction of variable sources in the bulge and ring regions 
given the data available. Kilgard et al. (2005) reported that 
29% of the sources show long-term variability based on the 
first two Chandra observations, consistent with our findings 
here. The fraction of variable sources in NGC 1291 is in gen- 
eral agreement with those in some typical early-type galax- 
ies, e.g., NGC 3379 and NGC 4278 (Brassington et al. 2008, 
2009). 

For the luminous nuclear point source and the three ULX 
candidates, we also examined their spectral variability by fit- 
ting the spectra in individual observations. The spectral vari- 
ability of the nuclear source is discussed in Section 4.2. For 
the ULX candidates, their spectra were well fit by an absorbed 
power-law model. Although all these ULX candidates are 
variable in flux, none of them shows significant (> 3a) vari- 
ability in its spectral parameters. The two non-transient ULX 
candidates (XIDS 26 and 55) still have ULX luminosities in 
individual observations. 
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Fig. 7. — Significance of flux variability vs. observed full-band luminos- 
ity for the 50 bulge sources (filled dots), 28 ring sources (open dots), and 3 
other sources (triangles). Dotted line indicates the > 3<r significance level of 
the flux variability. Note that we calculated flux variability for only sources 
detected in at least one individual observation; there arc insufficient data for 
variability study of the remaining sources. We studied spectral variability 
of the fev ULX candidates, which show no significant variability in spectral 
parameters (column density and power-law index). 

4.4. X-Ray Colors 

It is difficult to classify individual X-ray sources into differ- 
ent XRB types (e.g., LMXBs and HMXBs) given the limited 
X-ray data alone (e.g., Fabbiano 2006). However, we can use 
X-ray colors to separate the points sources into groups that are 
likely dominated by certain source types (e.g., Colbert et al. 
2004; Prestwich et al. 2003, 2009). We constructed an X- 
ray color-color plot following Prestwich et al. (2003). The 
X-ray colors are defined as HC = (H-M)/T (hard color) 
and SC = ( M-S)/T (soft color), where 5, A/, H , and T are 
the count rates in the 0.3-1. 0 keV, 1. 0-2.0 keV, 2.0-8. 0 keV, 
and 0.3-8. 0 keV bands, respectively. A Bayesian approach 
was adopted to calculate the colors and the associated er- 
rors following the method developed by Park et al. (2006). 
This approach provides a rigorous statistical treatment of the 
Poisson nature of the detected photons as well as the non- 
Gaussian nature of the error propagation. The color-color 
plot is shown in Figure 8. We consider a source to be sig- 
nificant in the plot (a black data point) if its average error of 
the soft and hard colors is smaller than the 3<x-clipped mean 
(0.18) of the errors for all the sources. The other sources 
(cyan data points) have relatively large color errors and are 
not included in the following analysis. A significant fraction 
(75% of the bulge sources and 65% of the ring sources) of 
the X-ray sources are located in the region that is likely dom- 
inated by LMXBs (Prestwich et al. 2003). The dominance of 
LMXBs is also suggested by the low SFR of NGC 1291. The 
SFR of NGC 1291 is ^ 20% of that of the Milky Way (e.g., 
Chomiuk & Povich 2011), and the Milky Way only hosts a 
few HMXBs with luminosities above our detection complete- 
ness limit (« 10 37 erg s“ l ; see Section 5.1) given its popu- 
lation studies (e.g., Grimm et al. 2002; Voss & Ajello 2010). 
Therefore, HMXBs should have a small contribution to the 
observed XRB population in NGC 1291 . 

The bulge and ring sources appear to have different distri- 
butions in the X-ray color-color space. A two-dimensional 
Kolmogorov-Smimov test 14 indicates that the probability of 
the two distributions being drawn from the same parent pop- 

14 We used IDL routine KS2d for this test; see 


FiG. 8. — X-ray color-color plot for the point sources. S, M, H, and T are 
the count rates in the 0.3-1 .0 keV, 1 .0-2.0 keV, 2.0-8.0 keV, and 0.3-8.0 keV 
bands, respectively. Black data points represent sources with relatively small 
color errors, while cyan data points represent sources with large errors and 
are less significant. The green and red tracks show the expected X-ray colors 
for absorbed power-law spectra with different power-law indices and col- 
umn densities. The blue ellipse indicates the area that is likely dominated by 
LMXBs (Prestwich et al. 2003). A significant fraction of the X-ray sources 
with small color errors are located in the region expected to be dominated by 
LMXBs. (A color version of this figure is available in the online journal.) 
ulation is 0.08. The 3<r-clipped mean soft and hard colors for 
the bulge population are 0.04 ± 0.02 and -0.14 ±0.02, and 
the mean values are 0. 1 1 ± 0.02 and -0. 1 5 ± 0.02 for the ring 
population. On average, the ring population appears to have 
a significantly harder soft X-ray color than the bulge popula- 
tion. The color-color plot is overlaid with a grid showing the 
expected locations of sources with different power-law spec- 
tra and absorption column densities, derived using the CXC’s 
Portable, Interactive, Multi-Mission Simulator (PIMMS). A 
generally harder soft color suggests that the ring sources are 
more obscured (larger Na) than the bulge sources, as would 
be expected in regions with ongoing/recent star formation. 

5. X-RAY LUMINOSITY FUNCTIONS 

We construct full-band XLFs for the bulge and ring regions 
using the merged Chandra observation. Although « 25-40% 
of the point sources are variable between the individual ob- 
servations, the XLFs derived from the merged observation 
should be stable against source variability based on previous 
studies of the XLFs for various systems (e.g., Grimm et al. 
2005; Zezas et al. 2007; Sell et al. 2011). 

5. 1 . Incompleteness Corrections and Background AGN 
Subtraction. 

To derive the XLFs correctly, it is essential to make proper 
corrections for the detection incompleteness and background 
AGN contamination. We performed simulations following 
Kim & Fabbiano (2004) to correct for the effects of incom- 
pleteness, including source detection limit and source confu- 
sion. Briefly, in each simulation we added a mock source at 
a random location of the event file of every observation using 
the MARX ray-tracing simulator, 15 created a merged event 
file and a merged full-band image, and then ran wavdetect 
to test if this additional source is detectable. The input X-ray 
luminosity of the source was randomly drawn from a power- 
law XLF with /? = 1 in a cumulative form, N(> Lx) = • 

http://www.astro.washington.edu/useiVyoachim/idVpy_idl.html. 

15 http://space.mit.edu/CXC/MARX/index.html 
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We assumed a power-law spectrum for the source with F = 
1 .7. The position of the source was randomly determined fol- 
lowing the r 1 / 4 law (de Vaucouleurs 1948). We note that the 
adopted luminosity and position distributions here do not af- 
fect the incompleteness corrections significantly, as we only 
aimed to derive the positional-dependent detection fractions 
at a given luminosity. We performed 60,000 simulations in 
total, and the incompleteness corrections were calculated as 
a function of X-ray luminosity for the bulge and ring re- 
gions, respectively, utilizing all the simulated sources in each 
region. We also excluded the central 10"-radius area from 
the bulge region when computing XLFs as the incomplete- 
ness corrections are not reliable in this crowded region. The 
resulting 90% (50%) completeness limits are ft 1.5 x 10 37 
(ft 7.1 x 10 36 ) and ft 2.2 x 10 37 (ft 8.5 x 10 36 ) erg s" 1 for 
the bulge and ring populations, respectively. The complete- 
ness limits in the bulge are lower because of the longer effec- 
tive exposure time in the bulge and the smaller background 
within the extraction regions of bulge sources. The XLFs 
were only calculated down to the 50% completeness lim- 
its, below which the Eddington bias (Eddington 1913) could 
have affected the measured X-ray luminosities. There are 64 
and 61 sources with luminosities above these completeness 
limits in the bulge and ring regions, respectively; ft 98% of 
the sources have more than 10 full-band counts and thus the 
Eddington bias should have minimal effect (e.g., Bauer et al. 
2004). Compared to the XLFs presented in Irwin et al. (2002) 
and Kilgard et al. (2005), which were based on the first two 
Chandra observations, we are probing the XLFs about twice 
as deep. 

To correct for background AGN contamination, we did not 
attempt to classify AGNs from the detected sources. Instead, 
we calculated the expected AGN luminosity function, which 
was then deducted from the observed XLFs. We computed the 
AGNs flux distribution using the Gilli et al. (2007) AGN PS 
model, and we normalized the AGN surface density to the ob- 
served value in the ft 4 Ms Chandra Deep Field-South, which 
is ft 10000 deg -2 at a flux limit of ft 3.2 x 10 -17 erg cm -2 s _1 
(Xue et al. 201 1). A 20% uncertainty was included for this es- 
timation that was attributed to cosmic variance (e.g., Luo et al. 
2008). There are 4.2 ± 0.8 and 21 .5 ± 4.3 AGNs expected 
in the bulge and ring regions, respectively. The ring XRB 
population is heavily contaminated (ft 35%) by background 
AGNs. To derive the differential XLFs of the point sources 
in NGC 1291, we calculated the expected number of back- 
ground AGNs in a given luminosity bin and then subtracted it 
from the incompleteness-corrected number of X-ray sources 
in this bin. The cumulative XLFs were constructed by inte- 
grating the differential XLFs. 

The incompleteness-corrected and background subtracted 
XLFs for the bulge and ring populations are displayed in Fig- 
ure 9. The 1 a errors were determined by combining the Pois- 
son errors of the number of sources in each luminosity bin and 
the 20% uncertainty of the AGN surface density. 

5.2. Parameterization of the X-Ray Luminosity Functions 

The shapes of XLFs for different XRB populations (mainly 
LMXBs and HMXBs) have been studied extensively in the 
literature (see e.g,, Fabbiano 2006 for a review). For LMXBs, 
the XLF in the bright end ( L x > a few 10 37 erg s _1 ) can 
usually be parameterized with broken power laws. The 
break luminosity is ft (2-5) x 10 38 erg s" 1 , which may 
be related to the Eddington luminosity of neutron star bi- 
naries (e.g., Sarazin et al. 2000; Finoguenov & Jones 2002; 


Kim & Fabbiano 2004). Below the break luminosity, the 
XLFs have differential power-law slopes of 1. 8-2.2 (e.g., 
Kim & Fabbiano 2004); above the break, the XLFs cut off 
sharply, with differential slopes of > 2.8. Note that the 
XLFs for young (< 5 Gyr) LMXBs samples may not present 
this high-luminosity break (Kim & Fabbiano 2010). In the 
faint end (L x < a few 10 37 erg s" 1 ), there are limited data 
available for the XLFs of LMXBs. The XLFs of nearby 
early-type galaxies and the bulges of nearby spirals (includ- 
ing the Milky Way) suggest that there is a second break at 
ft (2-5) x 10 37 erg s" 1 , below which the XLFs appear flat- 
ter with differential slopes of ft 0.5-1 (e.g., Gilfanov 2004; 
Voss & Gilfanov 2007; Kim et al. 2009; Zhang et al. 201 1). 16 
The origin of this low-luminosity break is not clear, which 
could be related to magnetic braking during binary stellar 
evolution (e.g., Postnov & Kuranov 2005; Fragos et al. 2008) 
or different types of donor stars in high- and low-luminosity 
LMXBs (Revnivtsev et al. 2011). For HMXBs, the XLFs 
can be described with a simple power law with a differen- 
tial slope of ft 1.6 over a large luminosity range (ft 5 x 
10 36_i 0 40 erg S -1. Griram et al 2003). 

We fit the XLFs of NGC 1291 to parameterize their shapes. 
A single power law of the form 



and a broken power law of the form 

dN_ ifL x <L b 

dL * Wtr ifix>i » (3) 

were used, where A 1 and ,42 are the numbers of sources at 
the reference luminosites. We used the sherpa spectral fit- 
ting tool in the CIAO package to fit the XLFs in their differ- 
ential forms. It has been shown in previous studies that this 
is an efficient method to derive the characteristics of XLFs 
and it can yield consistent results with a maximum likelihood 
approach (e.g., Zezas et al. 2007; Kim et al. 2009; Sell et al. 
2011). The XLFs were read into SHERPA as typical X-ray 
spectra with the number counts binned in small luminosity 
bins, 6 log(Lx) = 0.05. Given the small numbers of sources 
in the bulge and ring regions, we used the Cash statistic in 
SHERPA for the model fits. The incompleteness corrections 
were supplied as an ancillary response file (ARF), and the 
AGN contribution was supplied as a table model in addi- 
tion to the intrinsic power-law model. The best-fit parameters 
with lcr errors are reported in Table 4, and the best-fit bro- 
ken power-law models are shown in Figure 9. For the single- 
power model, both XLFs have a power slope of ft 1.5. For 
the broken power-law model, both XLFs show a break lumi- 
nosity of ft 5 x 10 37 erg s ~ l . Above the break luminosity, the 
power slope is ft 2; below that, the bulge XLF has a slope of 
ft 1 while the ring XLF appears flatter, with a slope of ft 0.4. 
Note that the errors for the broken power-law parameters are 
relatively large. 

The best-fit parameters were derived using the Cash statis- 
tic, which does not provide a goodness-of-fit parameter an- 
alytically. In order to compare the quality of the power- 
law and broken power-law fits, we adopted the likelihood 

16 It has been suggested that LMXBs in different environments (e.g., in the 
field or globular clusters) probably exhibit notable difference in their XLFs 
(e.g., Kim et al. 2009). Such detailed discussion is beyond the scope of the 
current paper. 
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Fig. 9. — Differential [panel (a)] and cumulative [panel (b)] XLFs for the bulge (red filled dots or curve) and ring (blue open dots or curve) populations. The 
XLFs have been corrected for detection incompleteness and AGN contamination. In panel (a), the dashed curves represent the expected XLFs of background 
AGNs for the corresponding regions. The ring XRB population is heavily contaminated by background AGNs, resulting in large error bars for the XLF data 
points. The solid curves are the best broken power- law fits to the corresponding XLFs. In panel (b), the red dotted curves represent the 1 <t errors for the ring 
XLF, and the grey region represents the la errors for the bulge XLF. The cumulative XLFs do not appear to monotonically decrease with luminosity at a few 
places dus to the uncertainties introduced by the background subtraction. The XLFs were computed down to the ~ 50% completeness limits. (A color version of 
this figure is available in the online journal.) 


ratio test used in Zezas et al. (2007). For each differential 
XLF (ring/bulge), we simulated 10,000 mock XLFs from the 
best-fit power-law model, using the same luminosity bins and 
adding Poisson fluctuations to the source counts. The de- 
tection incompleteness and background AGN contamination 
were also applied to the mock XLFs. We then fit the simulated 
XLFs with both the power-law and broken power- law models. 
The ratio of the Cash statistic values of the best power-law 
and broken power-law fits for the simulated XLFs were cal- 
culated and compared to the ratio for the observed XLF; the 
rate of obtaining a higher ratio in the simulations represents 
the probability of getting an impro\ ed broken power-law fit by 
chance. For the bulge XLF, we found that the broken power- 
law model provides a statistically improved fit at the 89% con- 
fidence level; for the ring XLF the broken power-law model 
provides an improved fit at the 94% confidence level. There- 
fore, the broken power-law model appears to provide better 
fits to the data for both the bulge and ring populations. The 
broken power-law shapes of the XLFs suggest that both the 
bulge and ring XRB populations are dominated by LMXBs. 
We are probing the low-luminosity break at ^ 5 x 10 37 erg s~ l , 
with a slope of « 2 at the brighter end and a flatter shape at the 
fainter end; there are too few luminous XRBs in NGC 1291 to 
constrain the high-luminosity break at ^ (2-5) x 10 38 erg s" 1 . 
This result is consistent with that inferred from the X-ray col- 
ors (see Section 4.4). Given the arguments above, we pre- 
fer the broken power-law model for both the bulge and ring 
XLFs, and the following discussions are based on the broken 
power-law parameters. 

6. THEORETICAL MODELING OF THE X-RAY BINARY 
POPULATIONS 

PS modeling has been proven a valuable tool in order to 
study the physical properties of the XRB populations ob- 
served with Chandra in nearby galaxies. Belczynski et al. 
(2004) constructed the first synthetic XRB populations for 
direct comparison with the Chandra observed XLF of NGC 
1569, a star forming dwarf irregular galaxy. Linden etal. 
(2009, 2010) developed models for the HMXB and Be XRBs 
of the SMC, studying the XLF and the spatial distribution of 
this population and investigating the effect of electron-capture 


supemovae of massive ONeMg stellar cores. Fragos et al. 
(2008, 2009) performed extensive PS simulations for the 
modeling of the two old elliptical galaxies NGC 3379 and 
NGC 4278. They studied the relative contribution to the ob- 
served XLF from sub-populations of LMXBs with different 
donor and accretor types, and the effects of the transient be- 
havior of LMXBs. All these studies have been pivotal in in- 
terpreting Chandra observations of nearby galaxies. 

In this paper, we utilized a new grid of PS models to in- 
terpret the XLF characteristics of NGC 1291. In these mod- 
els, we focused on LMXBs and HMXBs formed in the galac- 
tic field as products of the evolution of isolated primordial 
binaries (i.e. we did not take into account any dynamical 
effects). We performed our simulations with Star Track 
(Belczynski et al. 2002, 2008), an advanced PS code that has 
been tested and calibrated using detailed mass transfer star 
calculations and observations of binary populations, and in- 
corporates all the important physical processes of binary evo- 
lution. 

In the development of our models we incorporate our cur- 
rent knowledge about the characteristics of the stellar pop- 
ulation in NGC 1291. This galaxy has a stellar mass of 
« 1.5 x 10 11 M q and a low current SFR of « 0.4 M 0 yr -1 
(Kennicutt et al. 2003). Wc adopted the Noll etal. (2009) 
model for the star-formation history of NGC 1291 including 
two exponential starburst events (see Section 1 ). The metallic- 
ity of both populations is assumed to be solar. Due to the sim- 
ilar shapes of the bulge and ring cumulative XLFs and their 
relatively large uncertainties, it is not plausible to compare 
models to the individual XLFs. We also note that the best fit 
results by Noll et al. (2009) about the star-formation history 
of NGC 1291 have no spatial information. Hence, our com- 
parison between PS models and the observed XRB population 
is limited to the combined (bulge plus ring) XLF of the galaxy. 

We created a grid of 48 PS models varying parameters that 
arc known from previous studies (Fragos et al. 2008) to affect 
the shape and normalization of the resulting synthetic XLF. 
Namely, we varied the common envelope (CE) efficiency, the 
initial mass function, stellar wind strength, and finally we al- 
lowed for various outcomes of the CE phase, in particular the 
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possible CE inspirals with Hertzsprung gap donors that ter- 
minate binary evolution barring the subsequent XRB forma- 
tion (Belczynski et al. 2007). The latter parameter, which is 
known to affect the evolution of massive binaries which are 
the progenitors of double compact objects (Belczynski etal. 

2007) , is for the first time explored in the context of XRB 
populations. The various parameter values used are shown in 
table 5. We note that in our calculations we combined <*ce and 
A into cne CE parameter, where A is a measure of the central 
concentration of the donor and the envelope binding energy. 
In the rest of the text, whenever we mention the CE efficiency 
<*ce, we refer in practice to the product «ce x A, effectively 
treating A as a free parameter (see Belczynski et al. 2008, for 
details). 

The synthetic XLFs were constructed based on the method 
described in Fragosetal. (2008), while transient systems 
were modeled according to Fragos et al. (2009). Furthermore, 
we calculated the errors in the synthetic XLF due to Poisson 
low number statistics (high luminosity end) and uncertain- 
ties in the bolometric correction factors (low luminosity end) 
(T. Fragos et al. 2012, in prep.; P. Tzanavaris et al. 2012, in 
prep.). In order to assess which are the preferred values of the 
different model parameters, we calculated the likelihood of 
the observed XLF given a model. Our statistical analysis takes 
into account not only the shape of the XLF (as in Fragos et al. 

2008) but the absolute normalization too. Our maximum like- 
lihood model has a low CE efficiency (acE « 0.1) and the 
“standard” stellar wind prescription (Belczynski et al. 2010, 
and references therein). We also found that there is an appar- 
ent degeneracy between the slope of IMF and whether or not 
we restrict the outcome of a CE event to a merger when the 
donor star is in the Hertzsprung gap. A flatter IMF results in 
a larger black hole XRJB (BH-XRB) population, while on the 
other hand, assuming that a CE with a Hertzsprung gap donor 
leads always to a merger, which prevents many of the binaries 
from evolving into a BH-XRB. Hence, a BH-XRB population 
consistent with observations could either result from a steeper 
IMF (slope of -2.7) and CE prescription where we allow for 
all possible outcomes, even with a Hertzsprung gap donor, 
or it could result from a flatter IMF (slope of -2.35) and a 
CE model where a CE event with a Hertzsprung gap donor 
would always lead to a merger. In fact, the former case corre- 
sponds to our maximum likelihood model, while the latter to 
the model with the second highest likelihood. 

Figure 10 shows the synthetic XLF from our maximum 
likelihood model in comparison with the observed XLF, as 
well as the different types of XRBs contributing to the XLF 
at different ranges of X-ray luminosity. There are two ar- 
tificial “jumps” in the model XLF at « 10 37 erg s“ l and 
« 8 x 10 38 erg s"' 1 , instead of the observed smoother power- 
law shape. In our models, we identify transient and persistent 
sources by comparing the calculated mass-transfer rate of the 
binary to a critical mass-transfer rate, below which the ther- 
mal instability develops. Furthermore, we strictly limit the 
accretion rate to the Eddington limit. These limiting mass- 
transfer rates in our modeling are responsible for the apparent 
jumps in the synthetic XLF. However, in nature there is no 
sharp transition between thermally stable and unstable disks, 
nor a precise limit in the highest accretion rate possible. Ac- 
cretion onto a compact object is a non-linear and much more 
complex process, which in reality can results in a smoother 
luminosity distribution with no sharp transitions. Finally, by 
assuming a solar metallicity for the stellar population of NGC 
1291 in order to be consistent with Noll et al. (2009) and in 


the absence of any stellar metallicity measurements, we im- 
pose in practice a maximum BH mass of « 15 A/ 0 . If instead 
part of the stellar population had a lower metallicity (e.g. 30% 
solar) then the maximum BH mass would increase to 30 Mq 
or more (Belczynski et al. 2010), which could smooth out the 
jump at the very luminous end of the XLF. 

We found that the BH-XRBs that jpopulate the high-end of 
the XLF (above % 2 x 1 0 38 erg s , Eddington limit for a 
hydrogen accreting neutron star) are transient systems with 
evolved companion stars, while at lower luminosities the 
XRB population consists of a mixture of transient and per- 
sistent XRBs with a neutron star accretor and a giant or a 
main sequence (MS) donor star. These findings and the gen- 
eral shape of the synthetic XLF are in agreement with the PS 
study by Fragos et al. (2008) of the two old elliptical galax- 
ies NGC 3379 and NGC 4278. This leads us to the con- 
clusion that an old LMXB population is dominating the ob- 
served XLF. However, there is an evident overabundance of 
super-Eddington (Lx > 2 x 10 38 erg s" 1 , the Eddington limit 
for a hydrogen accreting 1 .4 M 0 neutron star) XRBs, com- 
pared to the Fragos et al. (2008) models, which can be at- 
tributed to the enrichment of the younger XRB population 
with black hole accretors. This younger population is likely 
originated from the ring star- forming region. Our finding that 
the younger stellar population of the ring region enhances 
the high X-ray luminosity, super- Eddington BH-XRB pop- 
ulation, is consistent with the recent observational study by 
Kim & Fabbiano (2010). Kim & Fabbiano (2010) compared 
the combined XLF of LMXBs detected in Chandra observa- 
tions of post-merger elliptical galaxies with that of typical old 
elliptical galaxies and found that elliptical galaxies with signs 
of recent star formation host a larger fraction of luminous 
LMXBs (L x > 5 x 10 38 erg s _l ) than old elliptical galaxies. 

7. DISCUSSION 

7. 1 . The Bulge LMXB Population 

Galactic bulges host relatively homogeneous samples of 
old stellar populations, and thus they are ideal targets for 
the research of LMXB XLFs. PS modeling of the XLFs has 
been used to constrain the formation, evolution, and physical 
properties of the LMXB population (e.g., Fragos et al. 2008, 

2009). Recently deep Chandra exposures of nearby ellipti- 
cal galaxies and spiral bulges have pushed the XLFs down 
to limiting luminosities of around a few 10 36 erg s _l ; some 
best studied examples include Centaurus A (Voss & Gilfanov 
2006), M31 (Voss & Gilfanov 2007), M81 (Sell et al. 201 1), 
and samples of galactic bulges in Gilfanov (2004) and 
Kim et al. (2009). The derived XLFs are remarkably simi- 
lar, with a low-luminosity break at « 5 x 10 37 erg s" and a 
high-luminosity break or cut-off at « 5 x 10 38 erg s“ l . The 
bulge of NGC 1291 is dominated by LMXBs based on the 
XLF shape and the X-ray colors, though relatively weak star 
formation could still exist given the observed filaments of Ha 
gas (See Section 1). The bulge XLF of NGC 1291 reaches a 
similar threshold luminosity of 7 x 10 36 erg s^ 1 , and its shape 
is statistically consistent with previous findings. We obtain a 
low-luminosity break at « 5 x 10 37 erg s -1 ; a high-luminosity 
break is not constrained but is highly probable as there is only 
one bulge source with Lx > 5 x 10 38 erg s' 1 (one of the ULX 
candidates). The power-law slope above the low-luminosity 
break is « 2, and is « 1 below the break, typical for LMXB 
XLFs. Despite of the complex morphology of the galaxy, the 
bulge of NGC 1291 still hosts a population of LMXBs that 
is very similar to those in typical elliptical galaxies and spiral 
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Fig. 10. — Left panel: Synthetic XLF (black solid curve) for the maximum likelihood model. The gray shaded area corresponds to the 1 -a error of the XLF 
due to Pcisson statistics and uncertainties in bolometric correction. For comparison the observed combined (corresponding to both the bulge and ring regions) 
XLF (red solid curve) and the associated 1 - o error (red dotted curves) are drawn. Center panel: Contribution of different sub-populations by separating the 
XRBs into groups of systems with different donor stellar types. Right panel: Contribution of different sub-populations by separating the XRBs into groups with 
different accretor types. In the center and left panel, solid colored lines denote persistent XRB sub-populations, while dashed colored lines denote transient ones. 
(A color version of this figure is available in the online journal.) 


bulges (in terms of their XLFs) down to the limiting luminos- 
ity, suggesting that all LMXB populations are compatible. 

7.2. Comparison of the Bulge and Ring Populations 

X-ray observations of the XRBs can provide valuable con- 
straints on the underlying stellar populations, independent of 
what were obtained from the other wavelengths. Given the 
X-ray observations of NGC 1291 , we consider that the ring of 
this galaxy hosts a younger stellar population than the bulge, 
based on the arguments given below. 

1. There appears to be an overdensity of X-ray sources 
in the ring of the galaxy. The radial profile of the X- 
ray source surface density follows the NUV surface 
brightness distribution in the ring instead of the AT-band 
surface brightness distribution, suggesting that the ring 
XRBs are likely associated with a younger stellar pop- 
ulation (see Section 4.1). 

2. The ring XRB population has in general a harder soft 
X-ray color than the bulge population, likely due to in- 
creased obscuration in the ring associated with ongo- 
ing/recent star-forming activities (see Section 4.4). 

3. The bulge and ring combined XLF can be well mod- 
eled with a synthetic XLF that is in general agreement 
with the PS modeling of typical old elliptical galaxies. 
However, there appears to be an enhancement of super- 
Eddington XRBs, which could be associated a possible 
younger stellar population in the ring region (see Sec- 
tion 6). 

4. The XLF is an efficient probe of the relative age of stel- 
lar populations. It has been reported that young ellipti- 
cal galaxies host a larger fraction of luminous LMXBs 
and the XLFs appear flatter than those of typical old 
elliptical galaxies (Kim & Fabbiano 2010). For NGC 
1291, the slope of the broken power-law fit at the faint 
end (Table 4) for the ring appears flatter than that for the 
bulge, albeit with a large uncertainty, suggesting that 
there are more luminous X-ray sources in the ring. A 
comparison of the numbers of sources in given lumi- 
nosity bins for the bulge and ring populations is shown 


in Figure 1 la. The numbers of sources have been cor- 
rected for the detection incompleteness and background 
AGN contamination. The fraction of relatively lumi- 
nous sources ( L x > 5.5 x 10 37 erg s" 1 )) is 30^% for 
the bulge, and it is 54^[^% for the ring; the errors are 
Poissonian. 

The LMXBs of the ring likely belong to two popula- 
tions, a relatively old population that ages similarly to 
the bulge population and a relatively young population 
that formed with the enhanced star formation. To iso- 
late this young population, we simply assume that the 
source number counts of the old population scale with 
the total Af-band luminosity, and then the old popula- 
tion in the ring can be derived from the bulge source 
number counts and the ratio between the AT-band lumi- 
nosities of the bulge and ring (a factor of 3.6). We sub- 
tract this old population from the ring number counts 
and the residual is displayed in Figure 1 lb. The frac- 
tion of relatively luminous sources is for this 

remaining population in the ring. The difference in the 
fraction of luminous sources between the bulge and the 
ring appears more significant, and it is likely that the 
remaining XRB population is indeed younger. 

7.3. XRB Population of Galactic Rings 

As the two best studied galactic rings in the X-ray, the 
rings of NGC 1291 and the Cartwheel galaxy differ signif- 
icantly. The XRB population in the ring of the Cartwheel 
galaxy is dominated by ULXs; the XLF has a power-law 
slope of « 1.6 albeit at high luminosities (Lx ^ 10 39 erg s -1 ; 
Wolter & Trinchieri 2004), consistent with being dominated 
by luminous young HMXBs. For NGC 1291, there are only 
two ULXs in the ring, and the ring XLF shape and X-ray 
source colors suggest a LMXB dominated system. The ring 
of the Cartwheel galaxy clearly has much stronger ongo- 
ing star formation than the ring of NGC 1291 (total SFR of 
0.4 M 0 yr" 1 for NGC 1291 compared to 67.5 Mq yr~ 2 for the 
Cartwheel galaxy). 

Such a discrepancy is probably caused by different forma- 
tion and evolution scenarios of the galactic rings. The ring of 
the Cartwheel galaxy is likely produced by a recent merger 
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Fig. 11. — Numbers of X-ray sources in given luminosity bins for the bulge, ring [panel (a)], and bulge-subtracted ring [panel (b); see text] populations. The 
numbers of sources have been corrected for the detection incompleteness and background AGN contamination. The luminosity bins were chosen to be the same 
for all the populations, and the limiting luminosity is as 8.5 x 10 36 erg s _1 , The luminosities for the bulge data points were shifted by a small amount for display 
purpose. For several cases with only error bars and no data points shown in the plot, the numbers of sources are smaller than 0.1. The vertical dashed lines mark 
the approximate breaks of the XLFs, and are used to separate relatively luminous and relatively faint sources. (A color version of this figure is available in the 
online journal.) 


event, and the expanding density wave forms the starburst 
ring. Collisional rings are generally short-lived with a lifetime 
of the order of a few hundred million years (see the review of 
Moiseev & Bizyaev 2009 and references therein). The ring of 
NGC 1291 is probably created by secular gas accumulation 
at the resonance orbit due to the gravity torques of the bars. 
Enhanced star formation appeared in the ring as gas was be- 
ing gathered. However, this could a be a less energetic and 
longer-term event than those appeared in collisional ring sys- 
tems, resulting in the relatively mild star formation observed 
today. Adopting the Noll et al. (2009) star-formation history 
model (Section 1) and assuming the ring hosts the young stel- 
lar population with an e-folding timescale of 50 Myr, the SFR 
of NGC 1291 would be ~ 20 A/ 0 yr~ ] 200 Myr ago, compa- 
rable to the current SFR of the Cartwheel galaxy. 

8. SUMMARY 

We presented Chandra observations of the ring galaxy 
NGC 1291 and studied the XRB populations in the ring and 
bulge regions of the galaxy. The main results are summarized 
below. 

1 . We analyzed the four Chandra ACIS-S3 observations 
of NGC 1291, with a total exposure of % 180 ks. We 
have detected 169 point sources in the 0.3-8. 0 keV 
band, 75 of which are in the bulge and 71 in the ring. 
The full-band luminosities of the sources range from 
2.7 x 10 36 to 2.3 x 10 39 erg s ~ l . 

2. We studied basic properties of the X-ray sources. There 
appears to be an overdensity of X-ray sources in the 
ring of the galaxy. The radial profile of the X-ray source 
surface density resembles the NUV surface brightness 
distribution in the ring, suggesting that a significant 
fraction of the ring XRBs are associated with a rel- 
atively young stellar population. About 40% of the 
sources in the bulge and 25% of the sources in the ring 
show >3 a long-term variability in the full-band count 
rate. 

3. We studied the spectra of the nuclear source, which ap- 
pears to be a low-luminosity AGN with moderate ob- 
scuration, and there is also hot gas present in the nuclear 


region that produces soft excesses in the spectra. The 
spectra of individual observations were well fit with an 
absorbed power-law plus thermal plasma model. 

4. In the X-ray color-color plot, a significant fraction of 
the bulge (« 75%) and ring (« 65%) sources with rel- 
atively small color errors are located in the region ex- 
pected for LMXBs. On average, the ring XRB popu- 
lation has a harder soft X-ray color, suggesting that the 
ring XRBs are more obscured, probably associated with 
ongoing/recent star-forming activities in the ring. 

5. We constructed full-band XLFs for the bulge and ring 
regions, taking into account the detection incomplete- 
ness and background AGN contamination. The re- 
sulting 90% completeness limits are « 1.5 x 10 37 and 
» 2.2 x 10 37 erg s~ l for the bulge and ring populations, 
respectively. The ring XRB population is heavily con- 
taminated (« 35%) by background AGNs. Both XLFs 
can be fit with a broken power-law model. The shapes 
are consistent with those expected for populations dom- 
inated by LMXBs. The ring XLF appears flatter than 
the bulge XLF, indicating that the ring hosts a larger 
fraction of luminous sources and is likely associated 
with a younger stellar population. 

6. We successfully modeled the bulge and ring combined 
XLF with a PS model, utilizing our current knowledge 
of the galaxy properties. The model is consistent with 
those for typical old elliptical galaxies, suggesting that 
an old LMXB population is dominating the observed 
XLF. The relative enrichment of luminous sources of 
NGC 1291 can be attributed a younger XRB population 
originated from the ring star-forming region. 

7. We consider that the ring of NGC 1291 hosts a younger 
stellar population than the bulge based on the X-ray 
studies presented in this paper. 
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TABLE 1 

Chandra OBSERVATIONS OF NGC 1291 


Obs. ID 
(1) 

Start Date 
(2) 

Exp (ks) 

(3) 

Cleaned Exp (ks) 

(4) 

pi 

(5) 

795 

2000 Jun 27 

39.7 

37.6 

J. A. Irwin 

2059 

2000 Dec 7 

37.0 

22.8 

A. H. Prestwich 

11272 

2010 May 4 

70.0 

69.1 

G. Fabbiano 

11273 

2010 May 8 

50.0 

49.1 

G. Fabbiano 


Note. — Column 1: Chandra observation identification number. Column 
3: Observation start data. Column 3: Nominal exposure time. Column 4: 
Exposure time after removing background flares. Column 5: Name of the 
Principal Investigator. 


TABLE 2 

Chandra POINT-SOURCE CATALOG 


X-ray Coordinates Counts 


XID 

(1) 

02000 

(2) 

^2000 

(3) 

Off-Axis 

(4) 

drad 

(5) 

FB 

(6) 

FB Low Err 

(V) 

FB Up Err 
(8) 

SB 

(9) 

SB Low Err 
(10) 

SB Up Err 
(11) 

1 ... 

03 16 47.17 

-41 08 15.2 

4.6 

6.2 

29.9 

5.9 

7.2 

19.2 

4.8 

6.2 

2 ... 

03 16 47.62 

-41 05 34.2 

4.2 

5.9 

144.8 

13.0 

14.2 

106.9 

11.3 

12.6 

3... 

03 16 48.09 

-41 04 10.1 

4.6 

6.2 

6.6 

2.8 

4.2 

10.8 

-1.0 

-1.0 


Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, aicminutes, and arcseconds. (This table is 
available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. The full 
table contains 34 columns of information on the 169 X-ray sources.) 


TABLE 3 

Best-Fit Parameters for the Nuclear Source 


Observation 

X 2 /dof 

r 

NR.! 

£-2 8 




(X 10 22 cm -2 ) 

(x 10 39 erg s' 1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 


795 

2059 

11272 

11273 


0.84 

0.93 

1.21 

U1 


$1 

21+H 


i 8*4:1 

!:Si 


9 c t+O .08 

■SSh 

1>Z/ -0.23 


Note. — The spectra were fit with an absorbed power law plus thermal plasma 
model (wabS 1*apec+wabs2*POW). Nm,i and T are die absorption column den- 
sity and photon index for the power-law component, respectively. The temperature 
(0. 15 keV) and column density = 7.6 X 10 21 cm' 2 ) of the thermal component 

were fixed at the values derived from die stacked spectrum. The observed 2. 0-8.0 
keV luminosities are dominated by the power-law component. 


TABLE 4 

Best-Fit Parameters for the X-Ray Luminosity Functions 


Region Single Power-law Broken Power-law 


Sp Or Al 5bp Qi Lfa (x 10 37 erg s" 1 ) A2 

Bulge 143 1.6^-j 36^ 138 l,l + g| 2.0*j-f 4.9^ 26+j 

Ring 159 1.42} 18 3 153 6A -u <Cf 


Note. — Sp and 5b P are the Cash fit statistic values for the power-law and broken 
power-lav/ fit, respectively. See Equations (2) and (3) for the definition of the other 
parameters. 
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TABLE 5 

Model Parameters for NGC1291 


Parameter 

Notation 

Value 

Star Formation History 


Noll et al. (2009) 

Metallicity 

Z 

0.02 

Total Stellar Mass 

Mm 

1.5 x 10“ MO 

Binary Fraction 

Fbin 

50% 

IMF power-law slope 


-2.7 or -2.35 

CE Efficiency 

a CE 

0.1, 0.2, 0.3, or 0.5 

Stellar wind strength 

??wind 

0.25, 1.0, or 2.0 

CE during HG 


Yes or No 

Magnetic Braking 


Ivanova & Taam (2003) 



